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Abstract: We consider two ultrarelativistic shock waves propagating and colliding in five- 
dimensional Anti-de-Sitter spacetime. By transforming to Rosen coordinates, we are able to 
find the form of the metric shortly after the collision. Using holographic renormalization, we 
calculate the energy-momentum tensor on the boundary of AdS space for early times after 
the collision. Via the gauge-gravity duality, this gives some insights on bulk dynamics of 
systems created by high energy scattering in strongly coupled gauge theories. We find that 
Bjorken boost-invariance is explicitely violated at early times and we obtain an estimate for 
the thermalization time in this simple system. 
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1. Introduction 

With the advent of high energy colliders with collision energies in the TeV range, progress 
in understanding the problem of ultrahigh energy particle scattering now involves knowledge 
about the dynamics of Quantum-Chromodynamics (QCD) above or close to the deconfinement 
scale. For instance, a new state of matter of deconfined quarks and gluons with fluid-like 
properties seems to have been created in Au+Au collisions at total collision energies of ~ 
40 TeV at the Relativistic Heavy- Ion Collider (RHIC) [1-4]. The upcoming Large Hadron 
Collider (LHC) will collide protons at ~ 14 TeV and lead nuclei at up to ~ 1100 TeV. 
Although not a collider, the Pierre-Auger observatory probes collision energies much beyond 
that, at up to ~ 10 8 TeV. 

QCD is a complicated theory to solve, especially when asking about real-time dynamics 
at energies close to the deconfinement scale. Traditionally one could only resort to weak 
coupling approaches, which due to asymptotic freedom are guaranteed to work well for bulk 
systems having extremely high energy densities (though not necessarily for those probed by 
RHIC and the LHC). In particular, phenomenological descriptions of RHIC data by applying 
hydrodynamic simulations with low viscosity [5] or none at all [6-8] seem to suggest that at 
these energy densities, the system is not weakly coupled [9,10]. 

Recently, the conjectured duality between strongly coupled gauge theories and gravity [11] 
has opened up a new window for studying strong coupling phenomena in a range of different 
gauge theories (although a dual description of QCD remains elusive to date). Concerning 
real-time dynamics, a lot of progress has been made in calculating transport coefficients for 
hydrodynamics in strong coupling, such as shear viscosity [12]. However, these calculations 
probe the response of a static medium at finite temperature, which - while important for 
near-equilibrium dynamics - do not give insights into the early, far from equilibrium stages 
of a high energy particle collision. 
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Prom the gauge theory point of view this earliest stage following the collision has to 
describe the transition of the system to an equilibrium state with a well-defined temperature 
(if the system lives long enough), which is referred to as thermalization. The gravity dual 
picture of thermalization is the formation of a black hole in the bulk, its Hawking temperature 
being identified with the temperature in the gauge theory on the boundary of the Anti-de- 
Sitter (AdS) space. Thus the problem of thermalization in strongly coupled gauge theories 
becomes related to the problem of black hole formation, which has been noted before [13] 
(see also [14,15] for a different proposal on thermalization involving black holes). 

As an aside, it should be pointed out that thermalization in a high energy particle collision 
is not guaranteed, and even if once achieved, not easy to maintain due to the rapid (initially 
one-dimensional) expansion of the system. 

In an inspiring work Janik and Peschanski [16] showed that an expanding thermal system 
corresponded to a gravity dual where the black hole was moving in the fifth dimension. 
Subsequent work on gauge-gravity duality in expanding systems continues to clarify the near- 
equilibrium late-time behavior of high energy particle collision duals [17-27]. 

Little is known about the dynamics in strongly coupled gauge theories shortly after a 
high energy collision, e.g. far from equilibrium. Notable exceptions are studies assuming 
independence from longitudinal dynamics [28], treating one space dimension instead of three 
[29] and a characterization of the dual horizon structure following a collision of two already 
deconfined plasmas [30]. We differ from [30] by using a simpler model for the incident states, 
which allows us to calculate the energy-momentum tensor analytically. Differences of our 
approach to Refs. [28, 29] will be discussed below. 

In this article we shall consider the problem of two colliding infinite sheets of matter 
in Af = 4 SYM in the large iV c and large 't Hooft coupling limit. Via the gauge-gravity 
duality we reinterpret the problem as the collision of two shock waves on the boundary of a 
five-dimensional AdS space. Solving for the five dimensional metric and using holographic 
renormalization, one thus is able to extract information about the real time dynamics of the 
energy momentum tensor of the gauge theory after the collision. 

Our work is organized as follows: in Section § we construct the line-element for two 
colliding shock-waves in Rosen coordinates. Section || deals with holographic renormalization: 
we construct there the energy momentum tensor perturbatively in proper time and check that 
it is covariantly conserved and traceless. We provide a physical interpretation of our results 
in Section [|, and put them into the perspective of the literature. 

2. Setup and Solution 

Following [16], we consider for J\f = SYM in the strong coupling, large N c regime, an energy- 
momentum tensor (EMT) of the single-shock-wave-form 

T++ = 0, T+_ = , T__ = m Six') , T xx = T yy = (2.1) 
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in light-cone coordinates x^ = This form serves as a toy model for a large particle 

moving nearly at the speed of light along the x + direction with transverse energy density 
jgjp = The attribute "large" refers to the perpendicular directions, xj_ = (x,y). In the 
gravity dual description, in Fefferman-Graham coordinates [31] this corresponds to the AdS§ 
line element with a single shock-wave 

2 _ -2dx + dx~ + fiiz 4 5(x')dx~ 2 + <ixj_ 2 + dz 2 
ds — k , 

z A 

which is an exact solution of the Einstein equations with negative cosmological constant 
(normalized here to A = —6) 

Rfiv - T^g^R - ^g^v = o . (2.3) 



In the following we shall consider the collision of two shock- waves of the form (|2,1| , |2.2D , 
corresponding to an EMT before the collision (t < 0) of 

T++ = n 2 5{x + ) , T+_ = , T__ = m S(x~), T xx = T yy = . (2.4) 

This problem was posed originally in Ref. [16]. The simple question we want to address is: 
what is the form of the EMT after the collision, and in particular in the forward light-cone 
(x ± > 0)? 

This setup of particle collisions in strongly coupled gauge theories mirrors closely that 
of Kovner, McLerran and Weigert [32] who treated collisions using classical Yang-Mills dy- 
namics. Starting from a charge current = 5^(5(x =t )pi j 2(x_L) and solving the classical 
Yang-Mills equations D^F^ U = J u they showed that the resulting gauge fields (and hence 
the EMT) were functions of the product x + x~ only, and therefore "boost-invariant" in the 
sense of Bjorken [33]. Quantum fluctuations are expected to break this boost-invariance [34], 
since even tiny fluctuations are unstable to exponential growth [35]. Our model fl2.4|) in some 
sense corresponds to the simple case pi,2(xj_) = ^1,2 = const., but is extendable to a situation 
where the p's are taken from the Color Glass Condensate framework [36,37]. We shall report 
on this interesting generalization in a subsequent work [38]. 

In order to calculate the energy momentum tensor in the forward light-cone for strongly 
coupled M = 4 SYM, let us first use the coordinate transformation 1 [39] 

x + = u + -m9(v) z 4 + 2fil v 9 2 (v) z 6 , x~=v, z = z + 2/ii v 9(v)2 3 , (2.5) 
where 9{v) is the Heaviside step function, to bring Eq. ( |2.2[ ) into the so-called Rosen form, 



2 _ ~2dudv + cbc ± 2 + [l + QpLiz 2 v 9{v)] 2 dz 2 
~ z 2 [l + 2fnz 2 v 9{v)} 2 ' ! ~' <> ! 



In some of the expressions below positive powers of the ^-function will appear. For positive and negative 
argument the square of the ^-function is equivalent to the ^-function, so the only complication arises if the 
argument vanishes. However, it turns out that all expressions of the type n (x) are multiplied by x m with 
positive m, and thus this complication and the associated ambiguity of defining 6 n (0) is of no relevance here. 
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This form is advantageous since it implies a metric that is continuous across the light-like 
hypersurface v = 0. Indeed, since one can do a similar transformation for the second shock 
wave, the precollision line element can be written as a simple superposition of the two (c.f. 
[39]), 

2 _ -2dudv + dx ± 2 + [l + 6/xiz 2 v 6{v) + 6[i 2 z 2 u 9(u)] 2 dz 2 
z 2 [1 + 2^i 5 2 v 9(v) + 2/i 2 z 2 u 9(u)} 2 
Since the metric has to be continuous and piece- wise differentiable in these coordinates, 
all corrections to the above line element after collision have to be proportional to uv 9{u) 9{v) 
(see appendix |A] for a more detailed discussion) . The calculation of these corrections is most 
conveniently done by introducing the coordinates of proper time f = y/2uv and space-time 
rapidity fj = In-. In these coordinates the hyper-surface spanned by u = v = becomes 
f = and the condition 9(u)9(y) ^ translates into f being real and positive. Introducing 
\x = y / 2/iijU2 and Y = \ In ^ , we then make the following ansatz for the line element after 
the collision 

2 ^ -[1 + K(f,fj,z) ] df 2 + [1 + L(f, fj, ~z)] f 2 df, 2 + [1 + H(f, fj, z)] dx ± 2 



z 2 [1 + 2z 2 /xfcosh(y — fj)] 



2 



+ [1 + M(f, fj, z)] [1 + 6z 2 fif cosh(y - fj)] 2 dz 2 ^ 
z 2 [1 + 2z 2 /ifcosh(y — fj)] 2 

where K,L,H,M are functions that vanish at f = and have to be determined by solving 
the Einstein equations (2.3). 

Determining K, L, H, M for all values of f may be possible with existing numerical meth- 
ods [38,40,41]. Finding full analytical solutions is much harder, so we limit ourselves to the 
restricted regime of early times f < 1, For this regime we use a power series ansatz in f 
for the functions K, L, H, M and determine the coefficients by solving the Einstein equations 
order by order in proper time. This is readily done with GRTensor [42]. One finds 

K(f, fj, ~z) = c^ 2 f 2 z 4 + 182 + 10c y f 3g6 cosh[Y -fj]- 4 z 2 

3 3 

838 + 160ci + c\ + 45c 2 4 _ 4 _ 8 154 + 2ci 4 _ 4 _ 8 . ,_ 5 

H fx t z n t z cosh[2(r — 7])] + 0(t ) 

9 3 

L (f , f} } ~ z ) = ~ 16 + C V 2 f 2 z 4 + 94 + 2 V f 3 i 6 cosh[Y - fj] - W 
3 3 3 

+c 2/ x W - 8Q6 + 1QC V W cosh[2(y - fj)] + 0(f 5 ) 
15 

H(f, fj, z) = -2 M W - 22 + 2Cl j »W cosh[y - fj] - ^±^/i W 

3 3 

_ 16 + 2^4-8 + 8_2ci 4 „ 8 cogh[2(y _ _ }] + 0( . 5) 

M(f, 77, 5) = 16/z 2 f 2 z 4 - 244 ~ 4C V W cosh[y - fj] + H + ^2f 4 z 2 

+ 1076 + 4c 1 /iW + 824-8^^ w cosh[2(y _ _ }] + 0(fB) (2 g) 
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where c±,C2 are freely choosable integration constants. Amusingly, all terms of 0(f 4 z 2 ) in 
( |2.9D cancel for the choice c\ = —5; however, at higher orders there is no similar freedom, 
and e.g. the O ) terms cannot be canceled. It is straightforward, but a little tedious 
to calculate K, L, H, M to arbitrary order in f . We have performed the calculation up to 
(including) 0(f 10 ), but believe little insight can be gained by spelling out this solution here. 
The constants ci,C2 reflect a residual gauge freedom, 

g^u = Cd a g^ + g m d v i a + g av d^ a . (2.10) 

To exhibit the property of c\ as parameter of residual gauge transformations we focus here on 
the leading order, i.e., we consider only terms up to (including) order f 2 (with the exception 
of gfjfj, which has to be considered up to order f ). With the generator 

2 ~3 ~4 

^ = C X HT Z Q otherwige ( 2 _ n ) 

6 

the gauge transformation acts on the metric as follows: 



A 9rr = -4 ci/i W + 0(f 3 ) , £g g m = ^ Cl ^ TV + Oif 5 ) (2.12) 

Z Z o 



All other components of the metric are either not influenced at all, or only at order f 3 . 
Comparison of fl2,12|) with fl2,9|) to order f 2 reveals that the terms generated by the residual 
gauge transformation ( 2,1C| ) with generator ( p. 11 ) are precisely the ci-dependent terms in 



( |2.9| ), This shows clearly that the freedom to choose c\ corresponds to a residual gauge 
freedom of the (partially) gauge fixed metric (^^), and thus we should expect that physical 
quantities, like the EMT, are independent from c\. We shall demonstrate this in the next 
Section. A similar analysis applies to higher orders in f, but the corresponding generator of 
residual gauge transformations is considerably more complicated than (2.11). 



3. Holographic Renormalization 

Having determined the solution to the metric for short times after the collision in the previous 
Section, we focus now on extracting information about the gauge theory EMT in this Section. 
Holographic renormalization [43] gives a simple prescription to obtain the EMT once the 
metric is in the Fefferman-Graham form, 

, 2 _ ^ Z^h^dM 

n=0 

where i collectively denote the coordinates on the AdS boundary [z = 0), gij is the metric on 
the boundary (assumed to be Minkowski) and Tj,- the gauge theory EMT. We can achieve to 
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bring Eq. (|2.8|) with (2.9) into the form (3.1) by the coordinate transformation 



f = r + X>(r,^ 4+2 " 

n=0 

oo 

n=0 

oo 



+ ^a n (r,r / )z 3 + 2 ", (3.2) 



z = z 

n=0 

where the z = transformation in f, has been chosen such that gij maintains the form 

gijdx l dx^ = —dr 2 + T 2 drj 2 + dx 2 ^ . 
One finds that the matching to lowest orders requires 

oq(t, r)) = -2/iT coshfy - nl - -it^tfV - fj, 3 T 7 cosh[y - n] + 0(r 8 ) 

6 

*o(t,t?) = -a T a (T,?7) 

eo(T,ry) = -— -^^ao(r, ry) (3.3) 

and in turn leads to the EMT 

T rr = fi 2 r 2 - 3^ 3 r 5 cosh[y - rj] + (107 + 90 cosh[2(y - r/)]) + C(r 10 ) 

45 

= -3^ 3 r 6 sinh[y - 77] + — ^ 4 r 9 sinh[2(y - rj)} + 0(r 10 ) 
T ~ 2 T m = -3/x 2 r 2 + 21//r 5 cosh[y - rj] - -//r 8 (107 + 150 cosh[2(Y -rj)]) + C(r 10 ) 



5_ 
2~4 r 

(3.4) 



T xx = T yy = 2/iV - 12juV cosh[y -7]] + — /i 4 r 8 (107 + 144cosh[2(y - r/)]) + C(r 10 ) . 



Note that the gauge theory EMT is independent from the residual gauge parameters c\ , c%, 

Moreover, it obeys Tjf = and V^T^" = 0, as it should, where we recall that is the 
covariant derivative with respect to the metric g^. 

4. Physics Interpretation 

Our main result, Eq. (|3.4|), gives the energy momentum tensor for short times after the 
collision of two sheets of matter in a strongly coupled gauge theory. A few remarks are in 
order: a non- vanishing off-diagonal element of Tj r in r, r\ coordinates means there is a flow 
of energy in the direction i, so the EMT is not in its local rest-frame. Locally, the EMT 
may always be brought into its rest-frame by a Lorentz boost (e.g. r\ — > rj + 0). However, 
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the form of Eq. (3.4) is such that the EMT may not be brought into its rest-frame globally 



(i.e. for all rj simultaneously), since it is explicitly dependent on rj. Put differently, the 



EMT (3.4) is not boost-invariant in the sense of Bjorken. This is a major difference to the 
result found when treating the gauge interaction as classical Yang- Mills fields (see discussion 
in the introduction). It may still be possible that boost-invariance is recovered (at least 
approximately) at late times, where our solution breaks down. The full solution obtained 
in a lower dimensional model [29], where boost-invariance (violated by construction at early 
times) is restored at late times, seems to suggest that this is the case. 

Nevertheless, it is interesting to discuss the time behavior of the EMT for one particular 
rest-frame, e.g. "central rapidity" 77 = Y. Then the individual diagonal components of the 
EMT have the interpretation of (local) energy density (T TT ), effective longitudinal pressure 
(T~ 2 T m ) and effective transverse pressure (T xx = T yy ). Keeping in mind that the result ( |3.4j) 
is valid only for r > since it does not contain the original discontinuities (^j4|) , it is interesting 
to note that all components of T^ u are very small initially and grow only proportional to r 2 
(although with negative sign in the case of the longitudinal pressure). We believe that this 
may be due to the simplicity of our model, and more specifically to the absence of transverse 
(xj_) dynamics in our ansatz for /i in fl2.4|). Along the same lines, Ref. [29] found that T^ u 
would vanish for all times had one started with (|2.4| ) and completely eliminated the transverse 
dimensions. 

At first glance, our result ( |3.4| ) contradicts that of Ref. [28], where it was found that 
the energy density should behave as a constant for < t < 1. However, Ref. [28] based 
their analysis on the assumption of boost-invariance, which is violated in our case, and we 
have neglected transverse dynamics which could change the behavior of lim r _i.o T TT . It may 
thus be possible to reconcile our results with that of Ref. [28]. On the other hand, Ref. [28] 
forbid solutions of rising energy density by invoking the positive energy condition T^uj^oj 11 > 
0, where uj^ is a time-like vector (see also [16]). While this criterion is certainly valid in 
classical gravity it is somewhat questionable why it should apply to the EMT of the boundary 
quantum field theory. After all, it is well-known that quantum fields cannot always and 
everywhere satisfy all energy conditions [44], and even the weakest form of energy conditions, 
the averaged null energy condition, can be violated for conformally coupled quantum fields in 
3+1 dimensions in any conformal quantum state [45]. It may be confusing to find quantum 
effects in a classical gravity calculation: however, it should be pointed out that these effects 
appear only at the boundary of AdS§, where the description should be dual to a strongly 
coupled quantum field theory [11]. Indeed, the five-dimensional EMT in our calculation 
trivially fulfills the positive energy condition, as it should for a classical gravity calculation. 

The breaking of boost-invariance of the EMT Eq. ( |3,4j ) is a direct consequence of the 
ansatz (|2.4| ) for the incoming Shockwaves. Namely, the relevant part of the precollision line 
element is 

ds 2 ~ z 2 [fi 1 6(x-)dx- 2 + LL 2 5{x + )dx +2 ] , (4.1) 
which - in terms of r, 77 coordinates - has an explicit 77 dependence. It is somewhat surprising, 
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1/3 



Figure 1: The early time behavior of the individual components of the EMT at rj = Y. 



though, that one can obtain a boost-invariant line element when formally replacing 5(x) — > 
d x 5(x), or similar structures 2 . In such a case, we could not find a coordinate transformation 
similar to Eq. (|2,5| ) in order to bring the line element into the Rosen form. It is conceivable 



that it would amount to replacing v9(v) and u6(u) in Eq. (2/7) by their derivatives, in which 
case the line element would be manifestly boost-invariant and the EMT start from a finite 
value. However, it is not evident how to interpret the physical meaning of a collision of two 
shock waves given by derivatives of delta functions. 

In Fig. @ we show a plot of the energy density and effective pressures in the local rest-frame 
for early times, which suggest that our small time expansion converges rapidly up to times 
Of pV 3 T~ 0.4. At this time, the effective longitudinal pressure is still negative and therefore 
the system is clearly not in equilibrium. At times > 0.7 the expansion in powers of 

r seems to break down, possibly signaling the onset of a transition to an equilibrated state, 
which from hydrodynamics is known to behave as T TT ~ t~ 4 / 3 . 



It would be interesting to study in detail the horizon structure of Eq. (|2.8|) . Short of doing 



2 Ref. [29] proposed 6(x)/x 2 to obtain boost-invariance. 
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this, we may hope to learn something about the position of the horizon by invoking cosmic 
censorship, e.g. requiring that all singularities of the metric fl2.8| ) are hidden behind horizons. 
To this end, it is instructive to consider the minimum value of z where H(f,fj,z) = — 1 at early 
times. Approximating H by its first term, we expect a singularity to appear at z 4 < tt4-^, 
which would have to be hidden by a horizon in order not to violate cosmic censorship. Note 
that this implies the horizon is moving towards the boundary at small times, in contrast to 
the situation at late times studied in Ref. [16]. Given that in our case the energy density is 
initially rising and that for static systems the temperature is inversely proportional to the 
distance of the black brane to the boundary, this is not too surprising. As a caveat, it should 
be pointed out that z 4 ~ is to be understood only as a very rough estimate of the 
singularity location since for these values of z essentially all terms of Eq. (|2.9| ) become of the 
same order, even at small times. 

Assuming that a singularity really does appear at , how long would it take until 

fl T 

the EMT at the boundary z = "knows" about the formation of the black hole (i.e. the 
thermalization of the system)? We again estimate the time by assuming the information gets 
transported over the distance z with the speed of light, so f ~ z. Hence we obtain a crude 
estimate of the thermalization time as a function of the transverse energy density fx, 

Ttherm ~ , (4.2) 

where the non-trivial dimensionless prefactor is 0(1) in our simple estimate (from Fig.|l], one 
can extract Ttherm > 0.4/i -1 / 3 ). Taken at face value, this would imply extremely small times 
r therm 1 fm/c for modern colliders like RHIC or the LHC. However, it seems that numerical 
simulations will be necessary to provide a detailed study of thermalization and extract the 
dimensionless prefactor in ( |4.2| ) reliably. Moreover, in a more realistic model than the one 
considered here, we expect the dynamics in the transverse coordinates xj_ to modify this 
thermalization time. We plan to study this in the near future [38]. 
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A. Distributions and Energy Conservation 



In order to address distributional issues we make the following global Ansatz for the line 
element in Rosen coordinates 

— 2dudvgi(u, v, z) + (u 2 dv 2 + v 2 du 2 ) g2(u, v, z) + g%(u, v, z)(fxj_ 2 + gi(u, v, z)dz 2 



ds z 



(A.l) 



For the precollision line element we can write 

gt(u,v,z) = g 3 (u,v,z) = fi(u,z)fi(v,z), g 2 (u,v,z) = 0, g 4 (u,v,z) = f 2 (u, z)f 2 (v, z), 
fl{u,z) = [l + 2jiz 2 u 6{u)Y 2 , f 2 (u,z) = [l + 6flz 2 u 6{u)] 2 [l + 2^~z 2 u 9(u)]~ 2 , 

where we have set \x\ = fi 2 = fi for simplicity. For this line element, the Einstein equations 
fl2.3|) do not contain any terms of the form 5(u), S(v ) which would be singular at u = 0, v = 0. 
Nevertheless, one might be worried that when making an ansatz for the line element for 
u > 0, v > this would introduce "spurious" singularities. To test for this, we choose the 
ansatz 

9i(u,v,z) = fi{u,z)fi(v,z) + 0(uv) 

g 2 (u,v,z) = e(u)e(v)(f 5 (u,z) + f 5 (v,z)) + 0(uv) 

g 3 (u,v,z) = f l (u,z)fi(v,z) +uv6(u)6(v) (fe(u,z) + fe(v,z)) + 0(u 2 v 2 ) 

g 4 (u, v, ~z) = f 2 (u, S)f 2 (v, z) + uv9(u)9(v) (f 7 (u, ~z) + f 7 (v, z)) + 0(u 2 v 2 ), (A.2) 

where fs(u, z), fe(u, z), fj{u, z) are required to be non-singular at u = 0, and we neglected 
terms of higher order that result in regular terms in fl2.3|) . Requiring that all terms of the 
form 5(u),5(v) cancel in the Einstein equations gives the condition 



(h(v,z) + / 5 (0,£))24/zW%) 



2 (l + 6j2z 2 ve(v)) 



f 7 (v,z)-f 7 (0,z) 



(1 + 2p,z 2 v6(v)) 

-2 [f 6 (v, z) + / 6 (0, z)} (1 + 6j2~z 2 v9(v)) 2 = 0, (A.3) 



and likewise for v <-> u. If Eq. ( |A.3j ) is fulfilled, then the equations ( |2.3| ) are regular at 
u = 0, v = and can be conveniently solved in r, r\- coordinates. By recasting our solution 



( |2.8D into the form (|A.2|), we find 

h(v,z) + f 5 (o,s) 1 

f 6 (v,z) + f 6 (0,z) 
fr(v,z) + f r (0,z) 



lim 



L-K 



(l + 2flz 2 v9(v)) 2 u^o 2uv 



1 



lim 

(1 + 2jlz 2 v9(v)) 2 u^o 

(1 + 6Rz 2 v9(v)) 2 ,. 
— -— hm 

(1 + 2jiz 2 v9(v)) 2 u^o 



lz A fx 2 



+ 



H 



(l + 2flz 2 v9(v)) uv 



8z 4 ft 2 



72z 4 [i 2 



+ 



M 



(l + 2flz 2 v9(v)) (1 + 6flz 2 v9(v)) uv 



.(A.4) 



Inserting the expressions (2J3) with cosh[n(y — r/)] = 2 ™ T ^n (u n fj, n + v n fi n ) and [it = 2a^/wv 
and reinstating appropriate 9 functions for every appearance of u, v, we obtain power series' 
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in v for the above f5,fe,f7, respectively. Using these, we have verified a posteriori that our 
solution obeys Eq. ( |A.3| ) order by order in v (and by symmetry also in u), and hence all 
5-functions in the Einstein equations cancel to that order. Since (|2.8| ) is a solution to the 
Einstein equations for arbitrarily small proper times, this implies that we know the metric 
foriiCl,i)C 1, including u = 0, v = 0. 

As a consequence of that, we are able to calculate the EMT for small (positive and 
negative) x^ by repeating the holographic renormalization procedure of section || for the 
metric in the form ( A.l ) to Brinkmann coordinates. One finds 

T__ = --d 2 _a (x + ,x-) + —9{x + )9(x~) lim (A.5) 
2 x~ z^o 2z 4 



which using 



5 + ci 



a (x + ,x~) = -2p(x + 9(x + ) +x~6{x~)) -— L p 2 4x +2 x- 2 9(x + )9(x-) + 0(x +3 ,x' 3 ) 

lim ~ = — (8 + Cl )fl 2 x + x-9(x + )9(x-) + 0(x +2 ,x~ 2 ) (A. 6) 

2^0 2z 4 3 



becomes 



T__ = pS(x~) -4x +2 fl 2 9(x + )9{x-) + 0(x +3 ,x-). (A.7) 



Note that in order for holographic renormalization to be consistent, all terms proportional to 
z 2 in L — K have to cancel. We have explicitly verified this up to 0(t 10 ). By symmetry we 
have 

T++ = j2d(x + ) - 4x- 2 p 2 9(x + )9(x~) + C(x~ 3 , x + ), (A.8) 



and from (3.4) we can glean 

T+_ = ^(T TT -T- 2 T m ) = 8x+x-jl 2 9{x + )9{x~) + 0{x' 2 ,x +2 ). (A.9) 

This can be used to calculate the energy density in the more familiar t, z coordinates. The 
transformation is straightforward and one finds 

T 00 (t, z) = - (T++ + 2T + - + T ) (A.10) 

~ -p (5(x + ) + 5(x~)) - 2p 2 9(x + )9(x-) {x +2 + x~ 2 - 4x+x~) 

~ -=ft (S(t + z) + 5(t - z)) + 2p 2 9(t 2 - z 2 ) (t 2 - 3z 2 ) . 
v2 

Taking the result (to this order of calculation) at face value, for fixed time t the energy density 
deposited in the forward lightcone is positive at the collision point z = (corresponding to 
mid-rapidity r] = 0), while it becomes negative for some z 2 < t 2 , especially close to the 
(positive) ^-functions at z 2 = t 2 . We believe this is allowed for a quantum field theory, as 
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discussed in section f|, but concede that the physical interpretation of this result demands 
further study. Calculating the total energy at any given time t gives 

/oo _ 
dzT 00 (t,z) = V2fl, (A.11) 
-oo 

up to the accuracy of the approximation. Therefore, the total system energy is conserved 
after the collision, which serves as yet another check on our result. 
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